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Abstract 

The time-dependent CP asymmetry in B® — > J/ipK + K~ decays is measured us- 
ing pp collision data at y/s = 7TeV, corresponding to an integrated luminosity of 
l.Ofb -1 , collected with the LHCb detector. The decay time distribution is charac- 
terised by the decay widths IY and Lh of the light and heavy mass eigenstates 
of the B®-B® system and by a CP-violating phase <fi s . In a sample of 27617 
B® — > J/tpK + K~ decays, where the dominant contribution comes from B® — > J/ij)(f> 
decays, these parameters are measured to be 

4> s = 0.07 ± 0.09 (stat) ± 0.01 (syst) rad, 

r s = (r L + r H )/2 = 0.663 ± 0.005 (stat) ± 0.006 (syst) ps" 1 , 
Ar s = r L -r H = O.lOO ± 0.016 (stat) ± 0.003 (syst) ps" 1 , 

corresponding to the single most precise determination of <f> s , AT S and T s . The 
result of performing a combined analysis with B® — > J/i)jit + it~ decays gives 

0.01 ± 0.07 (stat) ± 0.01 (syst) rad, 
0.661 ± 0.004 (stat) ± 0.006 (syst) ps" 1 , 
0.106 ± 0.011 (stat) ± 0.007 (syst) ps" 1 . 

All measurements are in agreement with the Standard Model predictions. 
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1 Introduction 



The interference between B® meson decay amplitudes to CP eigenstates J/ipX directly 
or via mixing gives rise to a measurable CP- violating phase <p s . In the Standard Model 
(SM), for b —7- ccs transitions and ignoring subleading penguin contributions, this phase 
is predicted to be —2(3 S , where (3 S = arg (— V ts V* b /V cs V* b ) and are elements of the CKM 
quark flavour mixing matrix [I] . The indirect determination via global fits to experimental 
data gives 2(3 S = 0.0364 ± 0.0016 rad This precise indirect determination within the 
SM makes the measurement of <p s interesting since new physics (NP) processes could 
modify the phase if new particles were to contribute to the B®-B® box diagrams j3j[I] 
shown in Fig. [TJ 

Direct measurements of S using B® — > J/ip4> and B® — > J7-07r + 7r~ decays have been 
reported previously. In the B® — > J/0 channel, the decay width difference of the light 
(L) and heavy (H) B® mass eigenstates, Ar s = IY — Th, and the average _B°-decay width, 
T s = (IY + Th)/2 are also measured. The measurements of <p s and Ar s are shown in 
Table IT] 

This paper extends previous LHCb measurements in the B° s — > J/ip(f) j5] and 
B Q S — > J/ipir + 7[~ [6] channels. In the previous analysis of B® — > J/ip4> decays, the invariant 
mass of the K + K~ system was limited to ±12 MeV/c 2 around the 0(1020) mass (7j, which 
selected predominately resonant P-wave — > K + K~ events, although a small S-wave 
K + K~ component was also present. In this analysis the K + K~ mass range is extended 
to ±30 MeV/c 2 and the notation B° s — > J/ipK + K~ is used to include explicitly both P- and 
S-wave decays |8). In both channels additional same-side flavour tagging information is 
used. The data were obtained from pp collisions collected by the LHCb experiment at a 
centre-of-mass energy of 7TeV during 2011, corresponding to an integrated luminosity of 
l.Ofb -1 . 

This paper is organised as follows. Section [2] presents the phenomenological aspects 
related to the measurement. Section[3]presents the LHCb detector. In Sect.|4]the selection 
of B® — > J/ipK + K~ candidates is described. Section [5] deals with decay time resolution, 
Sect. [6] with the decay time and angular acceptance effects and Sect. [7] with flavour 

Table 1: Results for S and Ar s from different experiments. The first uncertainty is statistical 
and the second is systematic (apart from the DO result, for which the uncertainties are com- 
bined). The CDF confidence level (CL) range quoted is that consistent with other experimental 
measurements of <j) s . 



Experiment Dataset [ fb 1 ] Ref . 4> s [ rad ] AR [ ps 1 ] 
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Figure 1: Feynman diagrams for B®-B® mixing, within the SM. 




tagging. The maximum likelihood fit is explained in Sect. [8] The results and systematic 
uncertainties for the B® — > J/tpK + K~ channel are given in Sections [9] and 10, the results 
for the B° s — > J j ^-n^ i\~ channel are given in Sect. 11 and finally the combined results are 
presented in Sect. [121 Charge conjugation is implied throughout the paper. 



2 Phenomenology 



The B® — > J/ipK + K decay proceeds predominantly via 5°-)- J/i>cj) with the (f) meson 
subsequently decaying to K + K~ . In this case there are two intermediate vector particles 
and the K + K~ pair is in a P-wave configuration. The final state is then a superposition 
of CP-even and CP-odd states depending upon the relative orbital angular momentum of 
the J/ip and the <fi. The phenomenological aspects of this process are described in many 
articles, e.g., Refs. 13j 14 . The main Feynman diagrams for B° s — > J/ipK + K~ decays 
are shown in Fig. [2 The effects induced by the sub-leading penguin contributions are 
discussed, e.g., in Ref. [15]. The same final state can also be produced with K + K~ pairs 
in an S-wave configuration 16 . This S-wave final state is CP-odd. The measurement 



of 4> s requires the CP-even and CP-odd components to be disentangled by analysing the 
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distribution of the reconstructed decay angles of the final-state particles. 

In contrast to Ref. jH], this analysis uses the decay angles defined in the helicity 
basis as this simplifies the angular description of the background and acceptance. The 
helicity angles are denoted by Q = (cos9k, cos 9^, tph) and their definition is shown in 
Fig. [3] The polar angle 9k (9^) is the angle between the K + momentum and the 
direction opposite to the B® momentum in the K + K~ (fi + fi~) centre-of-mass system. 
The azimuthal angle between the K + K~ and fi + fi~ decay planes is (fh- This angle is 
defined by a rotation from the K~ side of the K + K~ plane to the /i + side of the 
plane. The rotation is positive in the n + fi~ direction in the B° s rest frame. A definition 
of the angles in terms of the particle momenta is given in Appendix [Aj 

The decay can be decomposed into four time-dependent complex amplitudes, Ai(t). 
Three of these arise in the P-wave decay and correspond to the relative orientation of the 
linear polarisation vectors of the J/ip and mesons, where % G {0, ||, _L} and refers to the 
longitudinal, transverse-parallel and transverse-perpendicular orientations, respectively. 
The single K + K~ S-wave amplitude is denoted by A^{t). 

The distribution of the decay time and angles for a B° s meson produced at time t — 
is described by a sum of ten terms, corresponding to the four polarisation amplitudes 
and their interference terms. Each of these is given by the product of a time-dependent 
function and an angular function fl3] 



— >- oc fk(P) ■ (!) 



dt dQ 

fc=l 

The time-dependent functions hk{t) can be written as 

h k (t) = N k e~ Tat [a k cosh (|Ar a t) + b k sinh (§Aiy) 

+ Cfc cos(Am s t) + d k sin(Am s t)], (2) 

where Am s is the mass difference between the heavy and light B° s mass eigen- 
states. The expressions for the f k (Q) and the coefficients of Eq. [2] are given in Ta- 
ble [2] (XT , 18 . The coefficients N k are expressed in terms of the at t — 0, from 



now on denoted as Aj. The amplitudes are parameterised by |v4j|e i(Sl with the con- 
ventions 5q = and \Aq\ 2 + |A||| 2 + \A±\ 2 = 1. The S-wave fraction is defined as 
F s = \A S \ 2 /(\A \ 2 + |A,|| 2 + \A ± \ 2 + \A S \ 2 ) = \A S \ 2 /(\A S \ 2 + 1). 
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Table 2: Definition of angular and time-dependent functions. 
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C cos(<5q 


-S s ) 


S sin(5o — 


Ss) 


cos(<5q — 


Ss) 


D sin(<5o — 


Ss) 



For the coefficients a*., . . . ,dk, three CP violating observables are introduced 
1 - IA| 2 „ 23(A) _ 23?(A) 



C 



1 + IAI 



S 



1 + IAI 



1 + IAI 



(3) 



where the parameter A is defined below. These definitions for S and C correspond to 



those adopted by HFAG 
symbol A^ T used in Ref. 



19 



19 



and the sign of D is chosen such that it is equivalent to the 
The CP-violating phase (fi s is defined by <fi s = — arg(A) and 



hence S and D can be written as 



S 



2|A 


sin <f) s 


1 + 


A 


2 



D 



2|A 


cos 




1 + 


A 


2 



(4) 



20 associated with each polarisation 



The parameter A describes CP violation in the interference between mixing and decay, 
and is derived from the CP-violating parameter 
state i 

q A 
V A 



A, 



(5) 



where (A\) is the amplitude for a P° (B®) meson to decay to final state i and the 
complex parameters p = {B s \Bl) and q = [B s \Bl) describe the relation between mass and 
flavour eigenstates. The polarisation states i have CP eigenvalue rji = +1 for i e {0, ||} 
and rji = — 1 for i G {_L, S}. Assuming that any possible CP violation in the decay is the 
same for all amplitudes, then the product rjiAi/Ai is independent of i. The polarisation- 
independent CP-violating parameter A is then defined such that Aj = r/jA. The differential 
decay rate for a B° s meson produced at time t = can be obtained by changing the sign 
of Cfc and dk and by including a relative factor \p/q\ 2 . 

The expressions are invariant under the transformation 



((f) s ,Ar s ,6 ,S\\,5±,S s ) i — > (tt - 4> s , -AT S , -S , -6\\,ir- S±, -6 s ) 
which gives rise to a two-fold ambiguity in the results. 



(6) 
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In the selected 7r + 7r invariant mass range the CP-odd fraction of B® — > J/ip7i + 7i' 



decays is greater than 97.7% at 95% confidence level (CL) as described in Ref. 21 1. As 
a consequence, no angular analysis of the decay products is required and the differential 
decay rate can be simplified to 



dT(B° s ->■ J/^tt+tT 



dt 

3 Detector 



ex h 7 {t). (7) 



The LHCb detector (22 is a single-arm forward spectrometer covering the pseudorapidity 
range 2 < r\ < 5, designed for the study of particles containing b or c quarks. The detec- 
tor includes a high precision tracking system consisting of a silicon-strip vertex detector 
surrounding the pp interaction region, a large-area silicon-strip detector located upstream 
of a dipole magnet with a bending power of about 4Tm, and three stations of silicon- 
strip detectors and straw drift-tubes placed downstream. The combined tracking system 
has momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, 
and impact parameter resolution of 20 |J.m for tracks with high transverse momentum. 



Charged hadrons are identified using two ring-imaging Cherenkov detectors 23 . Pho- 
ton, electron and hadron candidates are identified by a calorimeter system consisting of 
scintillating-pad and pre-shower detectors, an electromagnetic calorimeter and a hadronic 
calorimeter. Muons are identified by a system composed of alternating layers of iron and 
multiwire proportional chambers. The trigger consists of a hardware stage, based on in- 
formation from the calorimeter and muon systems, followed by a software stage which 



applies a full event reconstruction 24 



Simulated pp collisions are generated using Pythia 6.4 |25| with a specific LHCb 



configuration 261. Decays of hadronic particles are described by EvtGen 27 in which 



final state radiation is generated using Photos 28 . The interaction of the generated 
particles with the detector and its response are implemented using the Geant4 toolkit 29 
as described in Ref. 1 30 1 . 



4 Selection of B® — > J/ipK + K candidates 

The reconstruction of — > J/ipK + K~ candidates proceeds using the decays 
J/ip — > fi + fi~ combined with a pair of oppositely charged kaons. Events are first required 



to pass a hardware trigger 24 , which selects events containing muon or hadron candidates 



with high transverse momentum (pt). The subsequent software trigger 24 is composed 
of two stages, the first of which performs a partial event reconstruction. Two types of 
first-stage software trigger are employed. For the first type, events are required to have 
two well-identified oppositely-charged muons with invariant mass larger than 2.7GeV/c 2 . 
This trigger has an almost uniform acceptance as a function of decay time and will be 
referred to as unbiased. For the second type there must be at least one muon (one high-px 
track) with transverse momentum larger than 1 GeV/c (1.7GeV/c) and impact parameter 
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larger than 100 [xm with respect to the PV. This trigger introduces a non-trivial accep- 
tance as a function of decay time and will be referred to as biased. The second stage of 
the trigger performs a full event reconstruction and only retains events containing a 
pair with invariant mass within 120 MeV/c 2 of the Jfi/j mass [7] and which form a vertex 
that is significantly displaced from the PV, introducing another small decay time biasing 
effect. 

The final candidate selection is performed by applying kinematic and particle 
identification criteria to the final-state tracks. The Jfip meson candidates are formed 
from two oppositely-charged particles, originating from a common vertex, which have 
been identified as muons and which have pt larger than 500 MeV/ c. The invariant mass 
of the yU + /i~ pair, m(/i + yU~), must be in the range [3030, 3150] MeV/c 2 . During subsequent 
steps of the selection, m(/i + /i~) is constrained to the J/ip mass [7). 

The K + K~ candidates are formed from two oppositely-charged particles that have 
been identified as kaons and which originate from a common vertex. The K + K~ pair 
is required to have a p<r larger than lGeV/c. The invariant mass of the K + K~ pair, 
m(K + K~), must be in the range [990, 1050] MeV/c 2 . 

The B® candidates are reconstructed by combining the Jfif) candidate with the K + K~ 
pair, requiring their invariant mass m(J/il)K + K~) to be in the range [5200, 5550] MeV/c 2 . 
The decay time, t, of the B° s candidate is calculated from a vertex and kinematic fit that 
constrains the B® — > J/tpK + K~ candidate to originate from its associated PV |3l) . The 
X 2 of the fit (which has 7 degrees of freedom) is required to be less than 35. Multiple 
B® candidates are found in less than 1% of events; in these cases the candidate with 
the smallest \ 2 is chosen. B° s candidates are required to have decay time in the range 
[0.3, 14.0] ps; the lower bound on the decay time suppresses a large fraction of the prompt 
combinatorial background whilst having a negligible effect on the sensitivity to <p s . The 
kinematic fit evaluates an estimated decay time uncertainty, 07. Candidates with <j t larger 
than 0.12 ps are removed from the event sample. 

Figure [4] shows the m(J/i/jK + K^) distribution for events originating from both the 
unbiased and biased triggers, along with corresponding projection of an unbinned maxi- 
mum log-likelihood fit to the sample. The probability density function (PDF) used for the 
fit is composed of the sum of two Gaussian functions with a common mean and separate 
widths and an exponential function for the combinatorial background. In total, after the 
trigger and full offline selection requirements, there are 27617 ± 115 B® — > J/ipK + K~ 
signal events found by the fit. Of these, 23 502 ± 107 were selected by the unbiased trigger 
and 4115 ± 43 were exclusively selected by the biased trigger. The uncertainties quoted 
here come from propagating the uncertainty on the signal fraction evaluated by the fit. 

Figure^shows the invariant mass of the and K + K~ pairs satisfying the selection 



requirements. The background has been subtracted using the sPlot 32 technique with 
m(J/ipK + K~) as the discriminating variable. In both cases fits are also shown. For the 



di-muon system the fit model is a double Crystal Ball shape 33 . For the di-kaon system 



the total fit model is the sum of a relativistic P-wave Breit-Wigner distribution convolved 
with a Gaussian function to model the dominant (ft meson peak and a polynomial function 
to describe the small K + K~ S-wave component. 
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Figure 4: Invariant mass distribution of the selected B® — > J/ipK + K~ candidates. The mass of 
the n + n~ pair is constrained to the J/tp mass (7j. Curves for the fitted contributions from signal 
(dotted red), background (dotted green) and their combination (solid blue) are overlaid. 




Figure 5: Background subtracted invariant mass distributions of the (a) fjT~ ~ \i~ and (b) K + K~ 
systems in the selected sample of B® — > J/ipK + K~ candidates. The solid blue line represents 
the fit to the data points described in the text. 



5 Decay time resolution 

If the decay time resolution is not negligibly small compared to the B® meson oscillation 
period 27r/Am s w 350 fs, it affects the measurement of the oscillation amplitude, and 
thereby <p s . For a given decay time resolution, a t , the dilution of the amplitude can 
be expressed as V = exp(— of Am^/2) [34]. The relative systematic uncertainty on the 
dilution directly translates into a relative systematic uncertainty on <f> s . 

For each reconstructed candidate, at is estimated by the vertex fit with which the 
decay time is calculated. The signal distribution of a t is shown in Fig. [6] where the sPlot 
technique is used to subtract the background. To account for the fact that track parameter 
resolutions are not perfectly calibrated and that the resolution function is not Gaussian, 
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a, [ps] 



Figure 6: Decay time resolution, at, for selected B® — > J/tpK + K signal events. The curve 
shows a fit to the data of the sum of two gamma distributions with a common mean. 




Figure 7: Decay time distribution of prompt J/tp K + K~ candidates. The curve (solid blue) is the 
decay time model convolved with a Gaussian resolution model. The decay time model consists 
of a delta function for the prompt component and two exponential functions with different decay 
constants, which represent the B® — > J/ipK + K~ signal and long-lived background, respectively. 
The decay constants are determined from the fit. The same dataset is shown in both plots, on 
different scales. 



a triple Gaussian resolution model is constructed 

where d is a common small offset of a few fs, fj are event-independent resolution scale 
factors and ft is the fraction of each Gaussian component, normalised such that ^2 fi = 1- 
The scale factors are estimated from a sample of prompt fi + K + K~ combinations 
that pass the same selection criteria as the signal except for those that affect the decay 
time distribution. This sample consists primarily of prompt combinations that have a 
true decay time of zero. Consequently, the shape of the decay time distribution close to 
zero is representative of the resolution function itself. 



(t - df 



(8) 



S 



Prompt combinations for which the muon pair originates from a real J/ip meson have 
a better resolution than those with random muon pairs. Furthermore, fully simulated 
events confirm that the resolution evaluated using prompt J/if) — > decays with two 

random kaons is more representative for the resolution of B® signal decays than the purely 
combinatorial background. Consequently, in the data only J/ipK + K~ events are used to 
estimate the resolution function. These are isolated using the sPlot method to subtract 
the fi + fi~ combinatorial background. 

The background subtracted decay time distribution for J/ip K + K~ candidates is shown 
in Fig. [7] using linear and logarithmic scales. The distribution is characterised by a prompt 
peak and a tail due to J/i/j mesons from B decays. The resolution model parameters are 
determined by fitting the distribution with a decay time model that consists of a prompt 
peak and two exponential functions, convolved with the resolution model given in Eq. [8} 

The per-event resolution receives contributions both from the vertex resolution and 
from the momentum resolution. The latter contribution is proportional to the decay time 
and cannot be calibrated with the prompt J/ijjK + K~ control sample. When using a 
scale factor for the resolution there is an assumption that the vertex contribution and the 
momentum contribution have a common scale. This assumption is tested in simulations 
and a systematic uncertainty is assigned. 

The effective dilution of the resolution function is calculated by taking its Fourier 



transform calculated at frequency Am s 34 



/oo 
dt cos(Am s t) R(t;a t ). (9) 
-00 

Taking into account the distribution of the per-event resolution, the effective dilution for 
the calibrated resolution model is 0.72 ± 0.02. This dilution corresponds to an effective 
single Gaussian resolution of approximately 45 fs. The systematic uncertainty accounts 
for uncertainties due to the momentum resolution scale and other differences between the 
control sample and signal decays. It is derived from simulations. 

The sample used to extract the physics parameters of interest consists only of events 
with t > 0.3 ps. The observed decay time distribution of these events is not sensitive to 
details of the resolution function. Therefore, in order to simplify the fit procedure the 
resolution function for the final fit (described in Sect. [8]) is modelled with a single Gaus- 
sian distribution with a resolution scale factor, rt, chosen such that its effective dilution 
corresponds to that of the multiple Gaussian model. This scale factor is r t = 1.45 ± 0.06. 



6 Acceptance 

There are two distinct decay time acceptance effects that influence the B° s decay time 
distribution. First, there is a decrease in reconstruction efficiency for tracks with a 
large impact parameter with respect to the beam line. This effect is present both in 
the trigger and the offline reconstruction, and translates to a decrease in the B° s meson 
reconstruction efficiency as a function of its decay time. This decrease is parameterised 
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Figure 8: B® decay time trigger-acceptance functions obtained from data. The unbiased trigger 
category is shown on (a) an absolute scale and (b) the biased trigger category on an arbitrary 
scale. 



by a linear acceptance function Et(t) oc (1 + (3t), which multiplies the time dependent 
B° s J/il)K + K~ PDF described below. The parameterisation is determined using a con- 
trol sample of B ± — > J/%l)K ± events from data and simulated B° s — > J/ip(p events, leading 
to (5 = (—8.3 ± 4.0) x 10 _3 ps _1 . The uncertainty directly translates to a 4.0 x 10 _3 ps _1 
systematic uncertainty on F s . 

Secondly, a non-trivial decay time acceptance is introduced by the trigger selection. 
Binned functional descriptions of the acceptance for the unbiased and biased triggers are 
obtained from the data by exploiting the sample of B° s candidates that are also selected 
by a trigger that has no decay time bias, but was only used for a fraction of the recorded 
data. Figure [8] shows the corresponding acceptance functions that are included in the fit 
described in Sect. [Sj 

The acceptance as a function of the decay angles is not uniform due to the forward 
geometry of LHCb and the requirements placed upon the momenta of the final-state 
particles. The three-dimensional acceptance function, eq, is determined using simulated 
events which are subjected to the same trigger and selection criteria as the data. Figure [9] 
shows the angular efficiency as a function of each decay angle, integrated over the other 
angles. The relative acceptances vary by up to 20% peak-to-peak. The dominant effect 
in cos 9^ is due to the pp cuts applied to the muons. 

The acceptance is included in the unbinned maximum log-likelihood fitting procedure 
to signal weighted distributions (described in Sect. [I]). Since only a PDF to describe the 
signal is required, the acceptance function needs to be included only in the normalisation 
of the PDF through the ten integrals J dQeu(fl) fk{Q)- The acceptance factors for each 
event i, en(f^), appear only as a constant sum of logarithms and may be ignored in the 
likelihood maximisation. The ten integrals are determined from the fully simulated events 



using the procedure described in Ref. 35 
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Figure 9: Angular acceptance function evaluated with simulated B® — > J/ip(p events, scaled by 
the mean acceptance. The acceptance is shown as a function of (a) cos 8k, (b) cos^ and (c) 
<Ph, where in all cases the acceptance is integrated over the other two angles. The points are 
obtained by summing the inverse values of the underlying physics PDF for simulated events and 
the curves represent a polynomial parameterisation of the acceptance. 



7 Tagging the flavour at production 

Each reconstructed candidate is identified by flavour tagging algorithms as either a B® 
meson (q = +1) or a B° s meson (q = —1) at production. If the algorithms are unable to 
make a decision, the candidate is untagged (g = 0). 

The tagging decision, q, is based upon both opposite-side and same-side tagging al- 
gorithms. The opposite-side (OS) tagger relies on the pair production of b and b quarks 
and infers the flavour of the signal B° s meson from identification of the flavour of the 
other 6-hadron. The OS tagger uses the charge of the lepton (/i, e) from semileptonic b 
decays, the charge of the kaon from the b — > c — > s decay chain and the charge of the 
inclusive secondary vertex reconstructed from 6-hadron decay products. The same-side 
kaon (SSK) tagger exploits the hadronization process of the b(b) quark forming the signal 
B®(B®) meson. In events with a B° s candidate, the fragmentation of a b quark can lead to 
an extra s quark being available to form a hadron, often leading to a charged kaon. This 
kaon is correlated to the signal B° s in phase space and the sign of the charge identifies its 
initial flavour. 

The probability that the tagging determination is wrong (estimated wrong-tag proba- 
bility, rj) is based upon the output of a neural network trained on simulated events. It is 
subsequently calibrated with data in order to relate it to the true wrong-tag probability 
of the event, u>, as described below. 

The tagging decision and estimated wrong-tag probability are used event-by-event in 
order to maximise the tagging power, e tag T> 2 , which represents the effective reduction 
of the signal sample size due to imperfect tagging. In this expression e tag is the tagging 
efficiency, i.e., the fraction of events that are assigned a non-zero value of q, and T> — l — 2u 
is the dilution. 
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Table 3: Calibration parameters (po, pi,(rj) and Apo) corresponding to the OS and SSK taggers. 
The uncertainties are statistical and systematic, respectively, except for Apo where they have 
been added in quadrature. 



Calibration 


Pa 


Pi 


(v) 


Ap 


OS 


0.392 ±0.002 ±0.008 


1.000 ±0.020 ±0.012 


0.392 


0.011 ± 0.003 


SSK 


0.350 ±0.015 ±0.007 


1.000 ±0.160 ±0.020 


0.350 


-0.019 ±0.005 



7.1 Opposite side tagging 

The OS tagging algorithms and the procedure used to optimise and calibrate them are 



described in Ref. 36 . In this paper the same approach is used, updated to use the full 
2011 data set. 

Calibration of the estimated wrong-tag probability, rj, is performed using approxi- 
mately 250 000 B + — > J/ip K + events selected from data. The values of q and rj mea- 
sured by the OS taggers are compared to the known flavour, which is determined by the 



charge of the final state kaon. Figure 10 shows the average wrong tag probability in the 
B — > J/ipK ± control channel in bins of rj. For calibration purposes a linear relation is 
assumed 

Po + ~rr +Pi(v ~ (v)) , 

Apo (10) 
w(l)= Po ^ + Pi(.V - (v)) , 

where oj{rj) and oJ(rj) are the calibrated probabilities for wrong-tag assignment for B 
and B mesons, respectively. This parametrisation is chosen to minimise the correlation 
between the parameters po and p\. The resulting values of the calibration parameters po, 
Pi, Apo and (rj) (the mean value of 77 in the sample) are given in Table [3j The systematic 
uncertainties for p and p\ are determined by comparing the tagging performance for 
different decay channels, comparing different data taking periods and by modifying the 
assumptions of the fit model. The asymmetry parameter Apo is obtained by performing 
the calibration separately for B + and B~ decays. No significant difference of the tagging 



efficiency or of p x is measured (Ae tsig = (0.00 ± 0.10)%, Ap x = 0.06 ± 0.04). Figure [10 
shows the relation between u and rj for the full data sample. 

The overall effective OS tagging power for B® — > J/ipK + K~ candidates is 
£ta g £> 2 = (2.29 ± 0.06)%, with an efficiency of e tag = (33.00 ± 0.28)% and an effective 
average wrong-tag probability of (36.83 ± 0.15)% (statistical uncertainties only). 

7.2 Same side kaon tagging 

One of the improvements introduced in this analysis compared to Ref. [5] is the use 
of the SSK tagger. The SSK tagging algorithm was developed using large samples of 



simulated B® decays to D s 7r + and J/ip <p an d is documented in Ref. 37 . The algorithm 



preferentially selects kaons originating from the fragmentation of the signal B° s meson, and 
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Figure 10: Average measured wrong-tag probability (cj) versus estimated wrong-tag probability 
(77) calibrated on B + — > J/tf] K + signal events for the OS tagging combinations for the background 
subtracted events in the signal mass window. Points with errors are data, the red curve represents 
the result of the wrong-tag probability calibration, corresponding to the parameters of Table [3} 



rejects particles that originate either from the opposite-side B decay or the underlying 
event. For the optimisation, approximately 26 000 B° s ^r Djtt + data events are used. The 



same fit procedure employed to determine the B° s mixing frequency Am s 38 is used to 
maximise the effective tagging power e taug V 2 . 

The calibration was also performed using B° s — > D~tt + events and assuming the 



same linear relation given by Eq. [TO} The resulting values of the calibration parame- 
ters (po,pi, Apo) are given in the second row of Table [3j In contrast to the OS tagging 
case, it is more challenging to measure po and p\ separately for true B or B mesons at 
production using B° s — > D~n + events. Therefore, assuming that any tagging asymmetry 
is caused by the difference in interaction with matter of K + and K~ , Ap is estimated 
using B + — > J/ipK~, where the p and distributions of the OS tagged kaons are first 
reweighted to match those of SSK tagged kaons from a large sample of fully simulated 
B° s — > D~n + events. 

The effective SSK tagging power for B° s ->- J/^K + K~ events is e tSLg V 2 = (0.89±0.17)% 
and the tagging efficiency is e tllg = (10.26 ± 0.18)% (statistical uncertainties only). 



7.3 Combination of OS and SSK tagging 

Only a small fraction of tagged events are tagged by both the OS and the SSK algorithms. 
The algorithms are uncorrelated as they select mutually exclusive charged particles, either 
in terms of the impact parameter significance with respect to the PV, or in terms of the 
particle identification requirements. The two tagging results are combined taking into 
account both decisions and their corresponding estimate of 77. The combined estimated 
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Figure 11: Distributions of the estimated wrong-tag probability, tj, of the B® — > J/ipK + K~ 
signal events obtained using the sPlot method on the JfyK + K~ invariant mass distribution. 
Both the (a) OS-only and (b) SSK-only tagging categories are shown. 



wrong-tag probability and the corresponding uncertainties are obtained by combining the 
individual calibrations for the OS and SSK tagging and propagating their uncertainties 



according to the procedure defined in Ref. 36 . To simplify the fit implementation, 
the statistical and systematic uncertainties on the combined wrong-tag probability are 
assumed to be the same for all of these events. They are defined by the average values 
of the corresponding distributions computed event-by-event. The effective tagging power 
for these OS+SSK tagged events is St az V 2 = (0.51 ± 0.03)%, and the tagging efficiency is 
e tag = (3.90 ± 0.11)%. 

7.4 Overall tagging performance 

The overall effective tagging power obtained by combining all three categories is 
£tag£> 2 = (3.13 ± 0.12 ± 0.20)%, the tagging efficiency is e tag = (39.36 ± 0.32)% and the 



wrong-tag probability is to = 35.9%. Figure 11 shows the distributions of the estimated 



wrong-tag probability r] of the B° s — > J/ipK + K signal events obtained with the sPlot 
technique using m(J/ijj K + K~) as the discriminating variable. 



14 



8 Maximum likelihood fit procedure 



Each event is given a signal weight, Wi, using the sPlot [32] method with m(J/ip K + K~) 
as the discriminating variable. A weighted fit is then performed using a signal-only PDF, 
denoted by S, the details of which are described below. The joint negative log likelihood, 
C constructed as 

-ln£ = -a WilnS, (11) 

events i 

is minimised in the fit, where the factor a = J2i W»/ J2i is used to include the effect 
of the weights in the determination of the uncertainties |39|. 

8.1 The mass model used for weighting 

The signal mass distribution, S m (m(J/ip K + K~); m B o, a m , r 2 i, /i), is modelled by a double 
Gaussian function. The free parameters in the fit are the common mean, m#o, the width 
of the narrower Gaussian function, <x m , the ratio of the second to the first Gaussian width, 
r 2i , and the fraction of the first Gaussian, f\. 

The background mass distribution, B m (m(J/ijj K + K~)) is modelled by an exponential 
function. The full PDF is then constructed as 

V m = f s S m + (l- f s ) B m , (12) 

where f s is the signal fraction. Fig. [4] shows the result of fitting this model to the selected 
candidates. 

8.2 Dividing the data into bins of m(K + K~) 

The events selected for this analysis are within the m(K + K~) range [990, 1050] MeV/c 2 . 
The data are divided into six independent sets, where the boundaries are given in Table |4| 
Binning the data this way leads to an improvement in statistical precision by separating 
events with different signal fractions and the analysis becomes insensitive to correction 
factors which must be applied to each of the three S-wave interference terms in the 
differential decay rate (fs, /g, fio in Table [2|. These terms are required to account for 
an averaging effect resulting from the variation within each bin of the S-wave line-shape 
(assumed to be approximately uniform) relative to that of the P-wave (a relativistic Breit- 
Wigner function). In each bin, the correction factors are calculated by integrating the 
product of p with s* which appears in the interference terms between the P- and S-wave, 
where p and s are the normalised m(K + K~) lineshapes and * is the complex conjugation 
operator, 

/•m H 

/ ps* dm(K + K-) = C SP e- i0SP , (13) 

where [m L , m H ] denotes the boundaries of the m(K + K~) bin, C*sp is the correction factor 
and #sp is absorbed in the measurements of 8s — 5j_. The Csp correction factors are given 
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Table 4: Bins of m(K + K~) used in the analysis and the Csp correction factors for the S-wave 
interference term, assuming a uniform distribution of non-resonant K + K~ contribution and a 
non-relativistic Breit-Wigner shape for the decays via the <fi resonance. 



m(K + K~) bin [MeV/c 2 ] 


Csp 


990 - 1008 


0.966 


1008 - 1016 


0.956 


1016 - 1020 


0.926 


1020 - 1024 


0.926 


1024 - 1032 


0.956 


1032 - 1050 


0.966 



Table 5: Parameters of the common signal fit to the m(J/ipK + K ) distribution in data. 



Parameter 


Value 


m B o [MeV/c 2 ] 
a m S [MeV/c 2 ] 

h 

r 2 i 


5368.22 ±0.05 
6.08 ±0.13 
0.760 ±0.035 
2.07 ±0.09 



in Table |4| By using several bins these factors are close to one, whereas if only a single bin 
were used the correction would differ substantially from one. The effect of these factors 



on the fit results is very small and is discussed further in Sect. 10, where a different S- 
wave lineshape is considered. Binning the data in m(K + K~) allows a repetition of the 
procedure described in Ref. 40 to resolve the ambiguous solution described in Sect. [I] by 
inspecting the trend in the phase difference between the S- and P-wave components. 

The weights, Wi, are determined by performing a simultaneous fit to the m(J/if)K + K~) 
distribution in each of the m(K + K~) bins, using a common set of signal mass parameters 
and six independent background mass parameters. This fit is performed for m(J/if}K + K~) 
in the range [5200, 5550] MeV/c 2 and the results for the signal mass parameters are shown 
in Tabled 



8.3 The signal PDF 

The physics parameters of interest in this analysis are T s , AT S , \A \ 2 , \A±\ 2 , F$, 5\\, 6±, 
5s, 4> s , |A| and Am s , all of which are defined in Sect. [2j The signal PDF, S, is a function 
of the decay time, t, and angles, Q, and is conditional upon the estimated wrong-tag 
probability for the event, T], and the estimate of the decay time resolution for the event, 
a t - The data are separated into disjoint sets corresponding to each of the possible tagging 
decisions q G {—1, 0, ±1} and the unbiased and biased trigger samples. A separate signal 
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PDF, S q (t, fl\a t , rj; Z, N), is constructed for each event set, where Z represents the physics 
parameters and iV represents nuisance parameters described above. 

The Sg are constructed from the differential decay rates of B° s and B° s mesons described 



in Sect. 



where 



Denoting d ' r <^ff**> by X and d ' r( *£ff**> ^ *, then 

S = Sq 
q fsgdidn' 



(14) 



S-l 



so 



[ (1 - u) X(t, tt;Z)+cu X{t, Q,]Z)\® R(t; a t ) e t {t) e n (Q) 
[ u X(t, fi; Z) + (l- u) X(t, Q;Z)]<8> R(t\ a t )] e t (t) e n (fi) 
[ X(t, fl; Z) + X(t, n-Z)]® R(t; a t )] e t (t) e n (ft). 



(15) 



Asymmetries in the tagging efficiencies and relative magnitudes of the production rates for 
B° s and B° s mesons, as well as the factor \p/q\ 2 are not included in the model. Sensitivity 
to these effects is reduced by the use of separately normalised PDFs for each of the tagging 
decisions and any residual effect is shown to be negligible. 

All physics parameters are free in the fit apart from Am s , which is constrained to the 
value measured by LHCb of 17.63 ± 0.11 ps _1 |38|. The parameter <5 S — 5j_ is used in the 
minimisation instead of 5s as there is a large (90%) correlation between 5$ and 5±. 

In these expressions the terms u and u represent the wrong-tag probabilities for a 



respectively, and £1X6 cL function of 
as given in Eq. 10 The 



candidate produced as a genuine B® or B° s meson, 
rj and the (nuisance) calibration parameters (pi,po, (rj), Ap 
calibration parameters are given in Table [3] and are all included in the fit via Gaussian 
constraints with widths equal to their uncertainties. 

The expressions are convolved with the decay time resolution function, R(t; at) 
(Sect. [5]). The scale factor parameter, r t , is included in the fit with its value constrained 
by a Gaussian constraint with width equal to its uncertainty. The £t{t) and Sq(Q) terms 
are the decay time acceptance and decay-angle acceptance, respectively. The two different 
trigger samples have different decay time acceptance functions. These are described in 
Sect. E 

Since this weighted fit uses only a signal PDF there is no need to include the distri- 
butions of either the estimated wrong tag probability, rj, or the decay time resolution for 
each event, a t . The physics parameter estimation is then performed by a simultaneous fit 
to the weighted data in each of the m(K + K~) bins for each of the two trigger samples. 
All parameters are common, except for the S-wave fraction Fg and the phase difference 
— 5±, which are independent parameters for each range. 



9 Results for B° s -t J/^K + K~ decays 

The results of the fit for the principal physics parameters are given in Table [6] for the 
solution with Ar s > 0, showing both the statistical and the total systematic uncertainties 
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described in Sect. [TOl 

The statistical correlation matrix is shown in Table [7| The projections of the decay 
time and angular distributions are shown in Fig. 12 It was verified that the observed 



uncertainties are compatible with the expected sensitivities, by generating and fitting to 
a large number of simulated experiments. 

Figure 13 shows the 68%, 90% and 95% CL contours obtained from the two- 
dimensional profile likelihood ratio in the (Ar s , <f> s ) plane, corresponding to decreases 
in the log-likelihood of 1.15, 2.30 and 3.00 respectively. Only statistical uncertainties are 
included. The SM expectation 41 is shown. 

The results for the S-wave parameters are shown in Table [8] The likelihood profiles 
for these parameters are non-parabolic and are asymmetric. Therefore the 68% CL in- 
tervals obtained from the likelihood profiles, corresponding to a decrease of 0.5 in the 
log-likelihood, are reported. The variation of 5s — S± with m(K + K~) is shown in Fig. 14 
The decreasing trend confirms that expected for the physical solution with (p s close to 
as found in Ref. 



zero. 



40 



All results have been checked by splitting the dataset into sub-samples to compare 
different data taking periods, magnet polarities, £>°-tags and trigger categories. In all 
cases the results are consistent between the independent sub-samples. The measurements 
of 4> s , AT S and T s are the most precise to date. Both Ar s and <j) s agree well with the SM 
expectation [2,41 . 

These data also allow an independent measurement of Am s without constraining it 
to the value reported in Ref. [38). This is possible because there are several terms in 
the differential decay rate of Eq. [TJ principally /i 4 and h 6 , which contain sinusoidal terms 
in Am s t that are not multiplied by sin0 s . Figure 15 shows the likelihood profile as a 
function of Am s from a fit to the data where Am s is not constrained. The result of the 
fit gives 

Am s = 17.70 ± 0.10 (stat) ± 0.01 (syst) ps" 1 , 



which is consistent with other measurements 138 42 - 441 . 
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Table 6: Results of the maximum likelihood fit for the principal physics parameters. The first 
uncertainty is statistical and the second is systematic. The value of Am s was constrained to the 
measurement reported in Ref. [38]. The evaluation of the systematic uncertainties is described 
in Sect. M 



Parameter 


Value 


r s [ps- 1 ] 
Ar s [ps- 1 ] 

\A±\ 2 
\Ao\ 2 
5\\ [rad] 
S± [rad] 
4> s [rad] 
|A| 


0.663 ± 0.005 ± 0.006 
0.100 ±0.016 ±0.003 
0.249 ± 0.009 ± 0.006 
0.521 ±0.006 ±0.010 
3.30^ ±0.08 
3.07 ±0.22 ±0.07 
0.07 ±0.09 ±0.01 
0.94 ±0.03 ±0.02 



Table 7: Correlation matrix for the principal physics parameters. 





r s 


AF S 


\A±\ 2 


l^ol 2 




s± 


4> s 


|A| 




[ps- 1 ] 


[ps- 1 ] 






[rad] 


[rad] 


[rad] 




r s [ps- 1 ] 


1.00 


-0.39 


0.37 


-0.27 


-0.09 


-0.03 


0.06 


0.03 


Ar s [ps- 1 ] 




1.00 


-0.68 


0.63 


0.03 


0.04 


-0.04 


0.00 








1.00 


-0.58 


-0.28 


-0.09 


0.08 


-0.04 


\M 2 








1.00 


-0.02 


-0.00 


-0.05 


0.02 


5|| [rad] 










1.00 


0.32 


-0.03 


0.05 


5± [rad] 












1.00 


0.28 


0.00 


0s [rad] 














1.00 


0.04 


\M 
















1.00 
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Figure 12: Decay-time and helicity- angle distributions for B® — > J/tp K + K~ decays (data 
points) with the one-dimensional projections of the PDF at the maximal likelihood point. The 
solid blue line shows the total signal contribution, which is composed of CP-even (long-dashed 
red), CP-odd (short-dashed green) and S-wave (dotted-dashed purple) contributions. 
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Figure 13: Two-dimensional profile likelihood in the (Ar s , <j) s ) plane for the B® — > J/tpK + K 
dataset. Only the statistical uncertainty is included. The SM expectation of 
Ar s = 0.082 ± 0.021 ps" 1 and 4> s = -0.036 ± 0.002 rad is shown as the black point with er- 
ror bar (2,41 . 



Table 8: Results of the maximum likelihood fit for the S-wave parameters, with asymmetric sta- 
tistical and symmetric systematic uncertainties. The evaluation of the systematic uncertainties 
is described in Sect. Ell 



m(K + K~) 


bin [MeV/c 2 ] 


Parameter 


Value 


Cstat (asymmetric) 


C"syst 


990 


- 1008 


F s 




0.227 


+0.081,-0.073 


0.020 






5 S - 5± 


[rad] 


1.31 


+0.78,-0.49 


0.09 


1008 


- 1016 


Fs 




0.067 


+0.030,-0.027 


0.009 






5 S - 5± 


[rad] 


0.77 


+0.38,-0.23 


0.08 


1016 


- 1020 


Fs 




0.008 


+0.014,-0.007 


0.005 






5s - S± 


[rad] 


0.51 


+1.40,-0.30 


0.20 


1020 


- 1024 


Fs 




0.016 


+0.012,-0.009 


0.006 






5s - S± 


[rad] 


-0.51 


+0.21,-0.35 


0.15 


1024 


- 1032 


Fs 




0.055 


+0.027,-0.025 


0.008 






5s - <1l 


[rad] 


-0.46 


+0.18,-0.26 


0.05 


1032 


- 1050 


F s 




0.167 


+0.043, -0.042 


0.021 






5s - 5± 


[rad] 


-0.65 


+0.18,-0.22 


0.06 
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Figure 14: Variation of 5s — 5± with m(K + K~) where the uncertainties are the quadrature sum 
of the statistical and systematic uncertainties in each bin. The decreasing phase trend (blue 
circles) corresponds to the physical solution with <j) s close to zero and Ar s > 0. The ambiguous 
solution is also shown. 




Figure 15: Profile likelihood for Am s from a fit where Am s is unconstrained. 
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10 Systematic uncertainties for B° s — >• J/ij)K + K de- 
cays 

The parameters Am s , the tagging calibration parameters, and the event-by-event proper 
time scaling factor, r t , are all allowed to vary within their uncertainties in the fit. There- 
fore the systematic uncertainties from these sources are included in the statistical uncer- 
tainty on the physics parameters. The remaining systematic effects are discussed below 
and summarised in Tables [9j 10 and 11 



The parameters of the m(J/tjjK + K~) fit model are varied within their uncertainties 
and a new set of event weights are calculated. Repeating the full decay time and angular 
fit using the new weights gives negligible differences with respect to the results of the 
nominal fit. The assumption that m(J/ipK + K~) is independent of the decay time and 
angle variables is tested by re-evaluating the weights in bins of the decay time and angles. 
Repeating the full fit with the modified weights gives new estimates of the physics param- 
eter values in each bin. The total systematic uncertainty is computed from the square 
root of the sum of the individual variances, weighted by the number of signal events in 
each bin in cases where a significant difference is observed. 

Using simulated events, the only identified peaking background is from 
B° J/0 K* (892)° events where the pion from the i^*(892)° decay is misidentified as 
a kaon. The fraction of this contribution was estimated from the simulation to be at most 
1.5% for m(J/4)K + K-) in the range [5200,5550] MeV/c 2 . The effect of this background 
(which is not included in the PDF modelling) was estimated by embedding the simulated 
B° — > J/ip K* (892)° events in the signal sample and repeating the fit. The resulting vari- 
ations are taken as systematic uncertainties. The contribution of B° s mesons coming from 
the decay of B+ mesons is estimated to be negligible. 

Since the angular acceptance function, en, is determined from simulated events, it is 
important that the simulation gives a good description of the dependence of final-state 
particle efficiencies on their kinematic properties. Figure [16] shows significant discrepancies 
between simulated B° s — > J/ip<fi events and selected B® — > J/ipK + K~ data events where 
the background has been subtracted. To account for this difference the simulated events 
are re-weighted such that the kaon momentum distribution matches the data (re-weighting 
the muon momentum has negligible effect). A systematic uncertainty is estimated by 
determining en after this re-weighting and repeating the fit. The changes observed in 
physics parameters are taken as systematic uncertainties. A systematic uncertainty is 
included which arises from the limited size of the simulated data sample used to determine 
en- 

The lower decay time acceptance is included in the PDF using the binned functions 
described in Sect. |6j A systematic uncertainty is determined by repeating the fits with 
the bin values varied randomly within their statistical precision. The standard deviation 
of the distribution of central values obtained for each fit parameter is then assigned as 
the systematic uncertainty. The slope of the acceptance correction at large lifetimes is 
= (—8.3 ±4.0) x 10 -3 ps _1 . This leads to a 4.0 x 10 -3 ps _1 systematic uncertainty on r s . 
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Figure 16: Background-subtracted (a) kaon and (b) muon momentum distributions for 
Bg —7- J/tp K + K~ signal events in data compared to simulated B® — > J/ip<f) signal events. The 
distributions are normalised to the same area. A larger deviation is visible for kaons. 



The uncertainty on the LHCb length scale is estimated to be at most 0.020%, which 
translates directly in an uncertainty on T s and AT S of 0.020% with other parameters 
being unaffected. The momentum scale uncertainty is at most 0.022%. As it affects both 
the reconstructed momentum and mass of the B® meson, it cancels to a large extent and 
the resulting effect on T s and Ar s is negligible. 

The Csp factors (Table [1]) used in the nominal fit assume a non-resonant shape for 
the S-wave contribution. As a cross-check the factors are re-evaluated assuming a Flatte 
shape 45 and the fit is repeated. There is a negligible effect on all physics parameters 
except 5s~ S±. A small shift (approximately 10% of the statistical uncertainty) is observed 
in 5s — S± in each bin of m(K + K~), and is assigned as a systematic uncertainty. 

A possible bias of the fitting procedure is investigated by generating and fitting many 
simplified simulated experiments of equivalent size to the data sample. The resulting 
biases are small, and those which are not compatible with zero within three standard 
deviations are quoted as systematic uncertainties. 

A measurement of the asymmetry that results from CP violation in the interference 
between B Q S -B° S mixing and decay is potentially affected by CP violation in the mixing, 
direct CP violation in the decay, production asymmetry and tagging asymmetry. In the 
previous analysis [5] an explicit systematic uncertainty was included to account for this. 
In this analysis the fit parameter |A| is added, separate tagging calibrations are used for 
and B° s decisions, as well as separate normalisations of the PDF for each tagging decision. 
Any residual effects due to tagging efficiency asymmetry and production asymmetry are 
shown to be negligible through simulation studies. 

The measurement of Am s determined from these data alone without applying a con- 
straint has been reported in Sect. |9j The dominant sources of systematic uncertainty 
come from the knowledge of the LHCb length and momentum scales. No significant sys- 
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Table 9: Statistical and systematic uncertainties. 



uULIl 


1 s 


Ar 


\A , I 2 


I4J 2 
l^ol 


/in 


°_L 


rh 
<Ps 


IAI 




[P b J 


LP S \ 






[rad] 


ii acii 


ii aQi 




Stat, uncertainty 


0.0048 


0.016 


0.0086 


0.0061 


+0.13 
-0.21 


0.22 


0.091 


0.031 


Background subtraction 


0.0041 


0.002 




0.0031 


0.03 


0.02 


0.003 


0.003 


J/iPK*° background 




0.001 


0.0030 


0.0001 


0.01 


0.02 


0.004 


0.005 


Ang. acc. reweighting 


0.0007 




0.0052 


0.0091 


0.07 


0.05 


0.003 


0.020 


Ang. acc. statistical 


0.0002 




0.0020 


0.0010 


0.03 


0.04 


0.007 


0.006 


Lower decay time acc. model 


0.0023 


0.002 














Upper decay time acc. model 


0.0040 
















Length and mom. scales 


0.0002 
















Fit bias 






0.0010 












Quadratic sum of syst. 


0.0063 


0.003 


0.0064 


0.0097 


0.08 


0.07 


0.009 


0.022 


Total uncertainties 


0.0079 


0.016 


0.0107 


0.0114 


+0.15 
-0.23 


0.23 


0.091 


0.038 



Table 10: Statistical and systematic uncertainties for S-wave fractions in bins of m(K + K ). 



Source 


bin 1 


bin 2 


bin 3 


bin 4 


bin 5 


bin 6 




Fa 




F s 




F s 




Stat, uncertainty 


+0.081 
-0.073 


+0.030 
-0.027 


+0.014 
-0.007 


+0.012 
-0.009 


+0.027 
-0.025 


+0.043 
-0.042 


Background subtraction 


0.014 


0.003 


0.001 


0.002 


0.004 


0.006 


J/iJjK* background 


0.010 


0.006 


0.001 


0.001 


0.002 


0.018 


Angular acc. reweighting 


0.004 


0.006 


0.004 


0.005 


0.006 


0.007 


Angular acc. statistical 


0.003 


0.003 


0.002 


0.001 


0.003 


0.004 


Fit bias 


0.009 




0.002 


0.002 


0.001 


0.001 


Quadratic sum of syst. 


0.020 


0.009 


0.005 


0.006 


0.008 


0.021 


Total uncertainties 


+0.083 
-0.076 


+0.031 
-0.029 


+0.015 
-0.009 


+0.013 
-0.011 


+0.028 
-0.026 


+0.048 
-0.047 



tematic effect is observed after varying the decay time and angular acceptances and the 
decay time resolution. Adding all contributions in quadrature gives a total systematic 
uncertainty of ±0.01 ps _1 . 
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Table 11: Statistical and systematic uncertainties for S-wave phases in bins of m(K + K ). 



Source 


bin 1 


bin 2 


bin 3 


bin 4 


bin 5 


bin 6 




5s - Sj_ 


5s - £>_i_ 


5s - Sj_ 


5s - S± 


5s - Sj_ 


5s - Sj_ 




[rad] 


[rad] 


[rad] 


[rad] 


[rad] 


[rad] 


Stat, uncertainty 


+0.78 

-0.49 


+0.38 
-0.23 


+ 1.40 

-0.30 


+0.21 
-0.35 


+0.18 

-0.26 


+0.18 
-0.22 


Background subtraction 


0.03 


0.02 




0.03 


0.01 


0.01 


J/ipK*° background 


0.08 


0.04 


0.08 


0.01 


0.01 


0.05 


Angular acc. reweighting 


0.02 


0.03 


0.12 


0.13 


0.03 


0.01 


Angular acc. statistical 


0.033 


0.023 


0.067 


0.036 


0.019 


0.015 


Fit bias 


0.005 


0.043 


0.112 


0.049 


0.022 


0.016 


Csp factors 


0.007 


0.028 


0.049 


0.025 


0.021 


0.020 


Quadratic sum of syst. 
Total uncertainties 


0.09 

+0.79 

-0.50 


0.08 

+0.39 
-0.24 


0.20 

+ 1.41 

-0.36 


0.15 

+0.26 
-0.38 


0.05 

+0.19 

-0.26 


0.06 

+0.19 
-0.23 
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11 Results for — >• J/i/j7t + 7t decays 

The B® — > J/tfj7T + 7T~ analysis used in this paper is unchanged with respect to Ref. |6] 
except for: 

1. the inclusion of the same-side kaon tagger in the same manner as has already been 
described for the B® — > J /ipK + K~ sample. This increases the number of tagged 
signal candidates to 2146 OS-only, 497 SSK-only and 293 overlapped events com- 
pared to 2445 in Ref. [6]. The overall tagging efficiency is (39.5 ± 0.7)% and the 
tagging power increases from (2.43 ± 0.08 ± 0.26)% to (3.37 ± 0.12 ± 0.27)%; 

2. an updated decay time acceptance model. For this, the decay channel 
B° — y J/ipK* (892)° , which has a well known lifetime, is used to calibrate the de- 
cay time acceptance, and simulated events are used to determine a small relative 
correction between the acceptances for the B° — y J/ipK* (892)° and B® — y J/ipn + n^ 
decays; 

3. use of the updated values of T s and Ar s from the B® —y J/ipK + K~ analysis pre- 
sented in this paper as constraints in the fit for <f> s . 

The measurement of 4> s using only the B° s — y J/ifj7T + ii~ events is 

<p s = -0.14 ±° ;ll± 0.01 rad, 

where the systematic uncertainty is obtained in the same way as described in Ref. (6l. 
The decay time resolution in this channel is approximately 40 fs and its effect is included 
in the systematic uncertainty. 

In addition, the effective lifetime T^o^j^ n+n _ is measured by fitting a single exponen- 
tial function to the B® decay time distribution with no external constraints on T s and 
Ar s applied. The result is 

T B°->j/ipn+n- = 1-652 ± 0.024 (stat) ± 0.024 (syst) ps. 
This is equivalent to a decay width of 

nfo^JA^+TT- = °- 605 ± °- 009 ( stat ) ± °- 009 ( s y st ) PS" 1 , 

which, in the limit <p s = and |A| = 1, corresponds to Th- This result supersedes that 
reported in Ref. 46 . The uncertainty on the B° lifetime [7] used to calibrate the decay 
time acceptance is included in the statistical uncertainty. The remaining systematic un- 
certainty is evaluated by changing the background model and assigning half of the relative 
change between the fit results with and without the decay time acceptance correction in- 
cluded, leading to uncertainties of 0.011 ps and 0.021 ps, respectively. The total systematic 
uncertainty obtained by adding the two contributions in quadrature is 0.024 ps. 
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Table 12: Results of combined fit to the 5° -)• J/^K + K~ and 5° ->■ J/V'vr+vr- datasets. The 
first uncertainty is statistical and the second is systematic. 



Parameter 


Value 


r, [ps- 1 ] 
Ar s [ps- 1 ] 

IA)| 2 
S\\ [rad] 
S± [rad] 
4> s [rad] 
|A| 


0.661 ±0.004 ±0.006 
0.106 ±0.011 ±0.007 
0.246 ±0.007 ±0.006 
0.523 ±0.005 ±0.010 
3.32 toil ± 0.08 
3.04 ±0.20 ±0.07 
0.01 ±0.07 ±0.01 
0.93 ±0.03 ±0.02 



Table 13: Correlation matrix for statistical uncertainties on combined results. 
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[rad] 


[rad] 


[rad] 




r s [ps- 1 ] 


1.00 


0.10 


0.08 


0.03 


-0.08 


-0.04 


0.01 


0.00 


AT S [ps- 1 ] 




1.00 


-0.49 


0.47 


0.00 


0.00 


0.00 


-0.01 








1.00 


-0.40 


-0.37 


-0.14 


0.02 


-0.05 


\M 2 








1.00 


-0.05 


-0.03 


-0.01 


0.01 


<5|| [rad] 










1.00 


0.39 


-0.01 


0.13 


5± [rad] 












1.00 


0.21 


0.03 


(f>s [rad] 














1.00 


0.06 
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12 Combined results for B® — > J/i/jK + K and 
B® — ► J/iJj7r + 7r~ datasets 

This section presents the results from a simultaneous fit to both — > J/tpK + K~ and 
B° s J/^7r + 7r- datasets. The joint log-likelihood is minimised with the common param- 
eters being T s , AT S , <ft s , |A|, Am s and the tagging calibration parameters. The combined 



results are given in Table 12 The correlation matrix for the principal parameters is given 
in Table M 

For all parameters, except T s and Ar s , the same systematic uncertainties as presented 
for the stand-alone B° s — > J/tpK + K~ analysis are assigned. For Y s and Ar s additional 
systematic uncertainties of 0.001 ps -1 and 0.006 ps -1 respectively are included, due to 
the B® — > J/ip7T + 7i^ background model and decay time acceptance variations described 
above. 



2N 



13 Conclusion 



A sample of pp collisions at y/s = 7TeV, corresponding to an integrated luminosity of 
1.0 fb~\ collected with the LHCb detector is used to select 27 617 ± 115 B° s J/^K + R- 
events in a ±30 MeV/c 2 window around the 0(1020) meson mass (7J. The effective tagging 
efficiency from the opposite-side (same-side kaon) tagger is e e s = 2.29 ± 0.22% (0.89 ± 
0.18%). A combination of data and simulation based techniques are used to correct for 
detector efficiencies. These data have been analysed in six bins of m(K + K~), allowing the 
resolution of two symmetric solutions, leading to the single most precise measurements of 
4> s , T s and AT S 

<p s = 0.07 ± 0.09 (stat) ± 0.01 (syst) rad, 
T s = 0.663 ± 0.005 (stat) ± 0.006 (syst) ps -1 , 
Ar, = 0.100 ± 0.016 (stat) ± 0.003 (syst) ps -1 . 

The B® — > J/tpK + K~ events also allow an independent determination of 
Am, = 17.70 ± 0.10 ± 0.01 ps" 1 . 

The time-dependent CP-asymmetry measurement using B° s — > J/ipn+Ti^ events from 
Ref. |6) is updated to include same-side kaon tagger information. The result of performing 
a combined fit using both B® — > J/ipK + K~ and B® — > J/ip7r + 'K~ events gives 

S = 0.01 ± 0.07 (stat) ± 0.01 (syst) rad, 
T s = 0.661 ± 0.004 (stat) ± 0.006 (syst) ps -1 , 
Ar s = 0.106 ± 0.011 (stat) ± 0.007 (syst) ps -1 . 

The measurements of 4> s , AT S and T s are the most precise to date and are in agree- 
ment with SM predictions [2 41 . All measurements using B® — > J/ipK + K~ decays 



supersede our previous measurements reported in Ref. [5], and all measurements using 
B® — > J/?pn + n~ decays supersede our previous measurements reported in Ref. (6j. The 



B° s J/^7r + 7r _ effective lifetime measurement supersedes that reported in Ref. 46 . The 
combined results reported in Ref. |6] are superseded by those reported here. Since the 
combined results for T s and AT S include all lifetime information from both channels they 
should not be used in conjunction with the B s — > J/ipir + i{~ effective lifetime measure- 
ment. 
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A Definition of helicity decay angles 

The helicity angles can be defined in terms of the momenta of the decay particles. The 
momentum of particle a in the centre-of-mass system of S is denoted by pf . With this 
convention, unit vectors are defined along the helicity axis in the three centre-of-mass 
systems and the two unit normal vectors of the K + K~ and fi + n~ decay planes as 



-KK 



n KK 



Pk+ x Pk- 



'K+ 



X 



Pk- 



P/i+ + Pfi- 



n 



mi 



The helicity angles are defined in terms of these vectors as 
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